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P_| ■ Context. We present spectroscopic observations of 182 disk galaxies (96 in the cluster and 86 in the field environment) in the region 

of the Abell 901/902 multiple cluster system, which is located at a redshift of z ~ 0.165. We estimate dynamical parameters of the 
four subclusters and analyse the kinematics of spiral galaxies, searching for indications of ram-pressure stripping. Furthermore, we 
focus on dusty red galaxies as a possible intermediate stage in the transformation of field galaxies to lenticulars when falling into the 

O cluster. 
Aims. We obtained multi-object slit spectroscopy using the VLT instrument VIMOS. We carried out a redshift analysis, determined 
velocity dispersions using biweight statistics, and detected possible substructures with the Dressler-Shectman test. We exploited 
rotation curves from emission lines to analyse distortions in the gaseous disk of a galaxy, as well as HST/ACS images to quantify the 
morphological distortions of the stellar disk. 

Methods. The presence of substructures and non-Gaussian redshift distributions indicate that the cluster system is dynamically young 
and not in a virialised state. We find evidence of two important galaxy populations. Morphologically distorted galaxies are probably 
subject to increased tidal interactions. They show pronounced rotation-curve asymmetries at intermediate cluster-centric radii and low 
rest-frame peculiar velocities. Morphologically undistorted galaxies show the strongest rotation-curve asymmetries at high rest-frame 
velocities and low cluster-centric radii. Supposedly, this group is strongly affected by ram-pressure stripping due to interaction with 
the intra-cluster medium. Among the morphologically undistorted galaxies, dusty red galaxies have particularly strong rotation-curve 
asymmetries, suggesting that ram pressure is an important factor in these galaxies. Furthermore, dusty red galaxies have on average a 
bulge-to-total ratio that is higher by a factor of two than cluster blue-cloud and field galaxies. 

The fraction of kinematically distorted galaxies is 75% higher in the cluster than in the field environment. This difference mainly 
stems from morphologically undistorted galaxies, indicating a cluster-specific interaction process that only affects the gas kinematics 
but not the stellar morphology. Also the ratio between gas and stellar scale length is reduced for cluster galaxies compared to the field 
sample. Both findings could be explained best by ram-pressure effects. 

Results. Ram-pressure stripping seems to be an important interaction process in the multiple cluster system A90 1/902. Dusty red 
galaxies might be a crucial element in understanding the transformation of field disk galaxies into cluster lenticular galaxies. 
Conclusions. 
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1. Introduction ternal physical mechanisms responsible for these relations is 

. . . . challenging. Galaxies in cluster environments may differ from 

Ever since the pioneering work by| Dressler] (|1980|), growing ev- field galaxies ow ing to higher initial densities, leading to ev- 
idence has been found for multiple links between galaxy evolu- lier co ii apse . This could mean a difference in evolution due 
tion and environment. Early-type galaxies (elliptical and SO or to their << nature ". But also the influence of "nurture" can be 
lenticular galaxies) are predominant in dense regions like clus- important . Both Galaxy-galaxy interactions, such as tidal in- 
ter cores, whereas late-type galaxies are more abundant in less teract ions, major/minor mergers, or harassment (iMoore et al.l 
dense field environments. Correlations have also been found be- |J996), an d galaxy-clus t er int e ractions, such as ram- pressure 



tween galaxy density and colour (e.g. |Blanton et al.| |2005|) or st[lppmg dGunn & Gottl (19721 iKronberger eTaTI l2008h due to 

star-formation rate (e.g.|Verdugo et al.|2008|>. Active galactic nu- me I CM (Intra-clu ster medium), halo truncat ion dBalogh et al. 

clei (AGN) seem to be more frequent in moderate-density envi- RoPC l l B ekki et al.ll200lh. or ga s compression dBvrd & Valtonen 

ronments ( GilmouretalJ|2007 ; KauffmannetalJ|2003). i QQffclTonnesen & Brvan 2009i can nlav crucial roles in the evo- 



1990; Tonnesen & Brvan 2009) can play crucial roles in the evo 
Disentangling the relative importance of the internal and ex- lution of ga i axies . To decouple the different effects of environ- 

ment on galaxy evolution, it is necessary to analyse a variety 
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rates (SFRs) over a wide range of environments. 
Some properties of high-density environments remain puzzling. 
One of these is the large increase in the fraction of lenticular 
galaxies in local galaxy clusters compared to higher redshifts, 
accompanied by a proportional decrease in the spiral fraction 
dDressler et al.lfl997b . A possible scenario would be that spiral 
galaxies falling from the field into a cluster's gravitational po- 
tential are transformed into SO galaxies. Ram-pressure stripping 
of the infalling galaxies' gas due to the I CM or starvation due to 
the cluster potential (e.g. Balogh et al. 2000) might be important 
factors in this transformation and responsible for the truncation 
of star formation. Besides this disk fading there would need to 
be some mechanism that causes the bulges in spirals to grow, 
since bulges in SOs are more luminous than those i n spirals (e.g. 



Christlein & Zabludo"f3 l2004). Galaxy harassment dMoore et al.l 
19961) or tidal interactions could be such mechanisms and hence 



might contribute to the transformation process. The dominant 
evolutionary path is still unknown. 

In th eir analysis of 10 inter mediate-redshift clusters (z ~ 0.4 - 
0.6), iPoggianti et all d 1999b find a suppression of star formation 
in cluster members. Furthermore they propose a candidate for 
an intermediate phase during the spiral-SO transformation. They 
identify a significant fraction of dusty starburst galaxies (strong 
Balmer absorption, modest [Oil] emission), probably falling 
in from the surrounding field. The descendants of these galax- 
ies could be poststarburst (k+a) galaxies and passive spirals. 
They suggest that the resulting passive disk populations migh t 
be the progenitors of SO galaxies (see also Jones eta~1 12000). 
In agreement with this hvpofhesis lKocevski et alj d201 1 ) find an 
obscured starburst population in cluster environments more fre- 
quently than in the field using Spitzer-24jum imaging. They pro- 
pose that, due to a mix of mergers and harassment, a higher frac- 
tion of cluster and group members are experiencing this nuclear 
starburst activity compared to field galaxies, and might be trans- 
formed into lenticulars with strong bulges. 
Another intermediate phase d uring this proces s could be the 
class of "dusty red spirals" dWolf et alJ 120031) . Having star- 
formati on rates four time s lower at fixed mass compared to blue 
spirals dWolf et al.l l2009h and showing similar extinction lev- 
els, dusty red galaxies can be understood as the low specific- 
SFR tail of the blue cloud. They show only weak spiral struc- 
tures and appear predominantly in the stellar mass range of 
log(M,/M ) = [10, 11], where they represent over half of the 
star-f orming galaxies in the A901/902 cluster system ( Wol f et al.l 
2009). Such a rich dusty red population indicates that if star- 
formation quenching does indeed happen, it must be a slow pro- 
cess and that any morphological transformation is delayed. At 
log(M,/M ) < 10 dusty reds are rare, suggesting a fast transfor- 
mation in the low mass regime. 

The STAGES (S pace Telescope Ab ell 901/902 Galaxy Evolution 
Survey) project dGrav et aT1 l2C)09a) and its large amount of data 
(see Section 2) provides information including galaxy mor- 
phologies, SEDs, SFRs and stellar masses for a specific region 
on the sky - the multiple cluster system Abell 901/902. This 
system might still be in the process of formation and not yet 
be virialised. It is therefore an interesting laboratory to inves- 
tigate the interplay b etween galaxy evolution and environment. 
iGallazzi et al.l (2009) demonstrate the general suppression of star 
formation at higher galaxy number densities in the STAGES sys- 
tem, and find that a high fraction of the remaining star formation 
in the cluster is either obscured or dominated by old stellar pop- 
ulations. Therefore dusty star-forming galaxies are predominant 
at intermediate densities. But there are also exampl es of a lack of 
dependence of galaxy properties on environment: Maltbv ~et alJ 



d2012l) find for STAGES spirals no dependence of the stellar dis- 
tribution in the outer stellar disk on galaxy environment and con- 
clude that differences in this distribution are related to internal 
mechanisms or minor mergers. 

In this work we exploit spatially resolved spectra from the 
VLT (Very Large Telescope) instrument VIMOS (Visible 
MultiObject Spectrograph), co mplemented with a lready existing 
data from the STAGES survey dGrav et al.ll2009al) . The paper is 
organised as follows. Section 2 gives a short introduction into 
the STAGES project, reports on the target selection and obser- 
vation of our VLT- VIMOS spectroscopy and describes the data 
reduction process. Section 3 presents the redshift analysis and 
gives an overview of the subcluster properties. Section 4 repre- 
sents the main part and contains the kinematic analysis exploit- 
ing rotation curve asymmetries. The last subsection deals with 
the prominent role of dusty red galaxies. A summary is provided 
in Section 5. 

Throughout this paper, we assume a cosmology with Q,„ = 0.3, 
Q. A = 0.7 and H = TOkms-'Mpc -1 . 

2. Observations and data reduction 

2.1. The STAGES data set 

This work utilises data from the STAGES project. STAGES is a 
multi-wavelength survey that covers a wide range of galaxy lu- 
minosities and galaxy densities. 

The structure under scrutiny is the multiple cluster system 
Abell 901/902 at z ~ 0.165, which comprises 4 subclusters 
and has been the subject of V-band (F606W) Hubble Space 
Telescope/ Advanced Camera for Surveys (HST/ACS) imaging, 
covering a (f.5 x 0°5 (~ 5 x 5Mpc 2 ) area of the multi-cluster 
system. Additionally 17-band COMBO- 17 (Classifying Objects 
by Medium -Band Observations) , Spitzer 24 pm, XMM-Newton 
X-ray data dGilmour et al .112007 ) and gravitational lensing maps 
dHevmans et al. 20081) are available. Photometric redshifts with 
errors S z /(l + z) ~ 0.02 up to R = 24, spect ral energy 
distrib utions dWolf et al.ll2003l) and stellar masses dBorch et al.l 
2006) have been derived from the COMBO- 17 survey as well 
as star-formati on rates from the COMBO-17 UV and Spitzer 
24 fim data dBell et alj |2007). Data from the Galaxy Evolution 
Explorer (GALEX), 2 degree Field (2dF) spectrograph, and 
Giant Meterwave Radio Telescope (GMRT) are also available. 
An ov erview of the publicly available data is given in lGrav et al.l 
d2009al) . 

2.2. VLT observations 

The spectroscopic observations of the Abell cluster system 
A901/02 with the VLT instrument VIMOS in MOS (Multiobject 
spectroscopy) mode were completed between February 8 and 
March 10, 2010 (ESO-ID 384.A-0813, P.I. A. Bohm). Because 
of the large field of view, the optical path of VIMOS is split into 
four channels (also called quadrants). The field of view of these 
four channels is 7' x 8' each, separated by 2' gaps. The CCD area 
of one quadrant is 2048 (spatial axis) x 4096 (wavelength axis) 
pixels. The instrument provides an image scale of 07205/pixel. 
The grism of our choice was the high-resolution grism HR-blue 
to obtain spatially resolved spectra, including the [Oil] emis- 
sion line doublet at 3726/29 A. Depending on the slit position 
on the CCD, a different ~ 2050A-range of the spectral range 
([3700A,6740A]) is seen on the CCD. The HR-blue grism has 
an average dispersion of 0.5 1 A/pixel and a spectral resolution of 
R ~ 2000 at a slit width of 1'.'2. The seeing conditions covered 
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the range 0.53" < FWHM < 1.14" according to the Differential 
Image Motion Monitor (DIMM). 

Four VIMOS pointings (with ~ 1' overlap) were used to cover 
the 1/4 square degree A901/902 field across the infall regions to 
the cluster centres. Two exposures of 1800 seconds each were 
taken for each pointing. 

The spectroscopic targets from the STAGES catalogue were se- 
lected to have stellar mass M, > 1O 9 M , absolute magnitude 
Mb < -1 8, a star-forming SE P (selection criterion in the public 
catalogue iGrav et aD d2009bl) : sed_type > 2), photometric red- 
shift 0.1 < Zphot < 0.26, and a visually confirmed disk compo- 
nent on the ACS images with an inclination angle i > 30. This 
selection criterion resulted in ~ 320 cluster and ~ 160 field can- 
didates. 

Using MOS slits with tilt angles up to ±45°, aligned along the 
apparent major axis of the galaxies, a position angle interval 
of 90° was covered. The manual slit placement was performed 
using the ESO software VMMPS (VIMOS Mask Preparation 
Software). Eventually, 238 slits were successfully placed. Some 
objects were included in two or more set-ups due to lack of alter- 
nate targets within the given sky area. These galaxies will later 
be used for consistency checks of our data reduction and analy- 
sis. Hence, spectra of 215 different galaxies were obtained. 

2.3. Spectroscopic data reduction 

The basic reduction steps are performed using the ESO- VIMOS 
pipeline (version 2.2.3) together with Gasgano (version 2.3.0), 
a Java-based data file organiser developed and maintained by 
ESO. For additional reduction steps and further processing we 
use the software package MIDASQ 

First, we create a combined master-bias with the pipeline recipe 
vmbias. We apply the master-bias in the reduction of the flat 
field, arc lamp, and scientific exposures. We process the avail- 
able flat field exposures with the recipe vmspflat. This pro- 
duces a normalised master-screen-flat and a non-normalised flat 
field, called combined-screen-flat, which we use to determine 
the spectral curvature model. However, we do not utilise the 
master flat in further reduction steps, since "the HR-blue grism 
data show not yet understood reflections, particularly evident in 
the flat-fields and arc-lamp frames" (excerpt of ESO readme for 
MOS users). 

In a next step the recipe vmspcaldisp processes the arc lamp ex- 
posure and calculates the spectral distortion model, mentioned 
above. Additionally, a cleaning of bad pixels and cosmic rays 
is performed. An output table, called extraction table, contains 
the coefficients of the wavelength calibration. In most cases we 
could provide an extraction table generated from previous runs 
to the recipe as input, in order to iterate the modelling of the 
spectral distortions. The wavelength calibration spectra with slits 
in the upper third of the CCD result in a low number of arc lamp 
lines. We tackle this issue by adding Argon-lines to the default 
reference catalogue. 

To achieve a satisfying wavelength calibration even for slits with 
a large tilt angle, we change a few crucial parameters from their 
default values. We increase the size of the search window (key- 
word vimos. Parameters. extraction.window) around the expected 
arc line position to 30 pixels. We lower the threshold for detect- 
ing a spectral feature (vimos. Parameters. line. ident. level) in order 
not to miss weak arc lines. Furthermore, we change the order of 
the polynomial (vimos. Paramete rs.slit.orde r) for wavelength so- 



1 Midas, the Munich Image Data Analysis System is developed and 
maintained by the European Southern Observatory (ESO) 



lution modelling along the spatial axis within each slit from to 
1. 

Finally, the recipe vmmosobsstare reduces the science frames 
including cosmic cleaning, sky subtraction, wavelength calibra- 
tion, slit rectification and 2D-spectra extraction. 
An additional cosmic cleaning, a co-addition of the two ex- 
posures, and a lD-object-spectra extraction is realised subse- 
quently in MIDAS. 

2.4. Spectroscopic redshifts 

For 188 of 215 galaxies we could measure spectroscopic galaxy 
redshifts using emission lines. For another 12 galaxies we could 
determine the redshift using absorption lines. These galaxies, 
though selected as star forming, show no detectable emission 
lines in their spectra. For 15 objects (~ 7%) no redshift mea- 
surement is possible due to the signal-to-noise ratio being too 
low. Note that none of these are primary targets but are intstead 
"fill-up" objects. We show a redshift histogram of all 200 galax- 
ies in Figurefl] clearly peaking at the cluster centre at z ~ 0.165. 
The "bump" at z ~ 0.26 is probably a volume effect. 
Errors for galaxy redshifts can be estimated using the spectra 
of galaxies that have been observed in more than one mask 
(Milvang- Jensen et al.ll2008l) . For 16 galaxies we obtained two 
redshifts, for two objects three redshifts and for one object we 
even have four redshifts. For each galaxy and redshift, we cal- 
culate the difference between the redshift and the mean of the 
redshifts available for a certain object. To compensate for not 
calculating the differences with respect to the mean of the un- 
derlying distribution we multiply these differences by a scal- 
ing factor (> 1) inferred from Monte-Carlo simulations by 
Milvang-J ensen et alj d2008l) . 

We th en calculated the biweight estimator of scale dBeers et alj 
119901) of the 42 differences to obtain a robust estimator for 
the standard deviation of the error distribution. It yields 6 Z = 
0.00019 as the estimate of the typical redshift error. For the clus- 
ter redshift z ~ 0.165, this error corresponds to S z /(l + z) ~ 
0.00016 or47km/s in rest-frame velocity. 



120 



100 



80 



"Z. 60 



40 



20 








0.10 0.15 0.20 0.25 0.30 0.35 0.40 

Z 



Fig. 1. Redshift histogram of the 200 galaxies of our sample 
having determined redshifts, peaking at the cluster centre at 
z ~ 0.165. 
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3. Properties of subclusters 

3.1. Subcluster redshifts and velocity dispersions 

Measurements of velocity dispersion and mean redshift provide 
a way to estimate the mass of a cluster, and to decide about clus- 
ter membership for each galaxy. These calculations should be 
ma de using a r obust statistic. We here adopt the biweight statis- 
tic (iBeers et al.lll990h . 

Prev ious weak lensing an alysis of the Abell 901/902 cluster sys- 
tem ( Hevmans et al . 2008) indicated four subclusters in the dark 
matter distribution, namely A901a, A901b, A902 and the SW- 
group (see Figure We use the coordinates of the four dark 
matter peaks as subcluster centres in our following analysis. 
Furthermore, we apply the biweight estimators of location and 
scale to determine the mean redshift z sc i as well as the veloc- 
ity dispersion cr sc i of each subcluster. For datasets up to ~ 50 
galaxies the biweight sta tistic is superio r to other estimators and 
insensitive to outliers (Be ers et al.l ll990). 

We select member candidates for each subcluster as follows. 
First we apply a common filter in redshift space (0.15 < Zi < 
0.18), which corresponds to a ~ 4500km/s interval around the 
cluster redshift z so i. This already excludes most background and 
foreground field galaxies. We then assign each selected galaxy 
to the nearest subcluster, in terms of projected distance to the 
subcluster centre on the sky. 

For further analysis we transform the galaxy redshifts z, to re- 
cession velocities vt using the relativistic formula: 

u, ■ - c ■ ^5 wherec : speed of light. ( 1 ) 

(1+ZiY + l, 

The biweight location estimator of these velocities u, gives an 
initial approximation of the redshift of each subcluster, z sc i- To 
obtain the final rest-frame velocity dispersion tr sc i, we transform 
all velocities i>; of a particular subcluster into its rest-frame: 



With these rest-frame velocities i>, ml as input data, the velocity 
dispersion cr sc ] can be calculated as the biweight estimator of 
scale, whereas the biweight estimator of location refines the ini- 
tial guess of the subcluster redshift z$ch 

For all galaxies in the interval [-3cr sc ], +3<x sc i] the biweight 
statistic is iteratively calculated in a loop until convergence, i.e. 
until all input galaxies fall into the 3cr-interval. Accordingly, we 
define those galaxies to be cluster members of the respective sub- 
cluster. 

By generating bootstrap sampl es from the final set of galaxy ve- 
locities, y, resl (IBeers etaDll990b . we perform Monte Carlo simu- 
lations (2000 replications) to estimate the asymmetric 68% error 
bars. For each bootstrap sample we determine the biweight esti- 
mator of scale (velocity dispersion <x scl ) and of location (subclus- 
ter redshift z sc i). All these values combined simulate the underly- 
ing distributions. In a first approximation the mean and standard 
deviation of each distribution give confidence intervals. Since 
the distributions of the bootstrap scale estimators are usually not 
Gaussian, another biweight is performed. From the confidence 
intervals we infer the error bars of the biweight estimators of the 
original sample, by calculating the difference to the interval lim- 
its. 

We employ this procedure to determine the velocity dispersion 
and mean redshift of all four subclusters as well as their asym- 
metric error bars. The results are summarised in Table |T] The val- 
ues for the velocity dispersion agree with those derived by Gray 



Subcluster Member z scl cr scl 



galaxies N [km/s] 



A901a 


31 


1642+ 00006 


^ ul -203 


A901b 


34 


n ift98+ 0004 
U.lb28_ 00004 


1175 +ioo 


A902 


26 


1665+ 00006 

o.iuuj _ 00006 


7Q9+150 

' 7z, -203 


SW group 


15 


1689+ 00004 


552 +97 

JJ -158 



Table 1. Velocity dispersion and mean redshift of the 4 subclus- 
ters. Additionally the number N of 3cr-member galaxies of each 
subcluster is given. 



et. al (in prep.) from 2dF spectra. Figure |3]shows histograms of 
the peculiar velocities in subcluster restframe Uj M of the final 
member galaxies calculated for the given value of z sc i- Notably, 
the subcluster A901b has many objects in the wings of the distri- 
bution, which may really be part of another subcluster. We will 
discuss this matter quantitatively in Subsection 3.2. 




57mOO.OOs 30.00s 56m00.00s 9h55m30.00s 

RA (J2000) 



Fig. 2. The dark matter reconstruction of the multiple cluster sys- 
tem A901/902 in orange-red colour sc ale derived from weak- 
lensing analysis (Hevm ans et alJl2~008h . Cluster members as a 
result of the biweight analysis (Section 3.1) are indicated with 
black diamonds. The four big crosses indicate the subcluster 
centres of A901a, A901b, A902 and SW group (from top left 
to bottom right), respectivel.y The radii of the 4 big circles are 
equal to the virial radii rzoo H08 (see Table|2]). 



3.2. Substructure 

To quantify the existence of substr ucture in galaxy clusters many 
statistical tests have been devised. [Dressier & Shectman (1988) 
developed a method to estimate the significance and location of 
cluster substructure using galaxy redshifts and projected sky po- 
sitions. We implement this three-dimensional (velocity-position) 
test in our analysis and apply it to each subcluster. 
Starting from a list of (sub)cluster members with measured po- 
sitions (ra, dec) and velocities, we determine for each galaxy 
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Fig. 3. Histograms of peculiar velocities in the cluster rest-frame, 
for the 4 subclusters. 

the local mean velocity Di oca i and velocity dispersion cr loca i from 
the sample of N nn nearest neighbours on the sky. Following 
iPinknev et al.l d 19961) we used N nn = V^V for the number of near- 
est neighbours. This maximises the sensitivity to significant sub- 
structures while reducing the sensitivity to Poisson noise. Next, 
we compare these local values to the global mean velocity v and 
velocity dispersion <x of the whole subcluster (see Table [TJ. A 
measure for the deviation from the global values can then be de- 
fined for each galaxy as 



N m + l 



[(Blood - + (Clocal - 0") 2 J . 



(3) 



Eventually, we compute the cumulative deviation A = E<5, which 
serves as a statistic for quantifying the substructure. By com- 
paring the A statistic with a set of 1000 Monte Carlo simula- 
tions by randomly reshuffling the velocities of the subcluster 
members, we can quantify the statistical significance of a sub- 
structure by the fraction P of si mulations that result in a h igher 
A value than the observed one dDressler & Shectmanlfl988n . A 
small value of P corresponds to a high significance. The crite- 
rion P < 0.1 is commonl y used as an indicator for substructures 
(e.g. lPopesso et al.ll2007l) . 

We visualise the possible position of substructures in Figure |4] 
markers whose size is proportional to e 6 depict the distribution 
of member galaxies on the sky, thus quantifying the local devi- 
ation from the global kinematics. Additionally, the symbols are 
coloured according to the subcluster rest-frame peculiar velocity 
of a galaxy. Some large circles in a given area indicate a cor- 
related spatial and kinematic variation and may locate kinemati- 
cally cooler systems within the cluster (e.g. subgroups) as well as 
flows of infalling galaxies. It is noteworthy, that the A statistic is 
insensitive in situations where substructures are superimposed. 
It relies on some spatial displacement of the centroids. 
Except for the SW-group, each of the other 3 subclusters shows 
significant substructures. For A902 (lower left) a clump of galax- 
ies (shown in blue) indicates a substructure that might be associ- 
ated with the more blueshifted subcluster A901a. A901b shows 
the strongest subclustering, in terms of having a high probability 
for substructure in combination with the clearly non-Gaussian 
galaxy velocity distribution (see Figure |3). Nine member galax- 
ies of A901b have a delta-statistic 5 > 2. In comparison A901a 
has four, A902 has five and the SW group no such member 



galaxies. Accordingly, Hev mans et al.1 (120081) find for A901a a 
significant substructure in the dark matter distribution, associ- 
ated with the infalling X-ray group A901ff. 
The presence of substructures and non-Gaussian redshift distri- 
butions (see Figure [3} underline that the cluster system is dy- 
namically young and not in a virialised state, which would be 
characterised by a Gaussian galaxy veloci ty distribution (e.g. 
iMerrall & Henrikserj 120031: iNakamura 120001) and, as indicated 
by N-body si mulations, by a low mass fraction included in sub- 
structures (e.g. Sha w et al.l2006l) . The presence of such substruc- 
tures can lead to an overestimation of dynamical parameters, es- 
pecially when applying the virial theorem (see next subsection). 



3.3. Dynamical parameters of the subclusters 

Estimating dynamical parameters of a cluster is not an easy task 
and various methods are available. Applying the virial theorem 
to positions and velocities of cluster member galaxies is the stan- 
dard method for estimating the mass of a self-gravitating system. 
This approach assumes that the system under scrutiny is in dy- 
namical equilibrium and that the galaxies in a cluster trace the 
mass, which means that the number density is proportional to 
the mass density. The virial mass estimator is derived from the 
Jeans equation and relates the total kinetic energ y to the poten- 
tial en ergy (2T + U — 0). It is given by (see e.g. lCarlberg et al.l 
1996): 



M„ 



3cr 2 ttR dv <t : 



3D 



(4) 



where the three-dimensional velocity dispersion <tid and the 
three dimensional virial radius r„ can be calculated from observ- 
ables using the line-of-sight velocity dispersion <x (see Section 
6) and the projected virial radius R pv . For more details see 
ICarlberg et ail (1 19961) . 

One interpretation is that the radius of virialisation corresponds 
to the radius where the internal material is virialised, but the ex- 
ternal material is still plunging into the cluster. Some simulations 
indicate that this happ ens at r?no, th e radius where the mean inte- 
rior density is 200p c (White 2001). Accordingly, a widely used 
mass estimator is M200' the mass within r2oo- 
Note that the masses calculated using Equation|4]depend on the 
radial extent of the sample, which in turn is connected to the 
mean interior density p(r D ) inside the virial radius. The mean 
density p is scaled to the critical density p c (z) = 3i/ 2 (z)/87rG, 
the density for a flat universe at an epoch represented by redshift 
z. 

To determine these values we follow Carlberg e t al.l (1 19971) and 
first correct M„ for the surface pressure term (2T + U = 
is replaced by 2T + U - 3PV, accounting for not integrating 
the Jeans equation to infinity), which is the exte rnal pressure 
from matter outside the v irialised region (see e.g. iGirardi et al.l 
119981: iThe & White! 19861 for details'), by employing a typical re - 



npl< 

duction at r„ of 20% (ICastellano et al.ll201 it iRines et al.ll2003h . 
Then, by assuming a singular isothermal sphere model for the 
density -radius relation we extrapolate r v to r2oo: 



^200 = r„ 



P(r v ) 



200p c 



1/2 



(5) 



Subsequently, one can deter mine the correspond ing mass M2oo- 
Using N-body simulations, Perea et al ] (119901) show that the 
virial mass estimator is superior to other two-mass estimators 
(median mass estimator, projected mass estimator) and is less 
sensitive to anisotropies or subclustering. However, it is affected 
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Fig. 4. Results from the Dressler-Shectman (DS) test. The plots show the spatial location of the cluster member galaxies. The radii of 
the plotted circles are proportional to e s , where 6 is the DS measurement of the local deviation from the global velocity dispersion 
and mean recessional velocity, i.e. larger symbols correspond to a higher significance of residing in a substructure. The colours 
indicate the rest-frame velocity relative to the subcluster centre, ranging from red (receding; 1500km/s) to blue (approaching; - 
1500km/s). The black crosses indicate the subcluster centres. P is the probability of there being no substructure in the dataset. 
Except for the SW-group, each of the other 3 subclusters shows a high significance for substructure. The significant substructures 
and non-Gaussian redshift distributions indicate that the cluster system is dynamically young and not in a virialised state. 



by the presence of interlopers, i.e. unbound and non-virialised 
galaxies, or the existence of a mass spectrum. Each factor can in- 
troduce an overestimation of the cluster mass by factors as large 
as two to four. Thus, the calculated masses have to be considered 
as a upper limits. 

Another bias might be introduced by the presence of a velocity- 
morphology segregation. For many clusters, the star-forming 
late-type population with blue colours shows a larger velocity 
dispersion than the red and passive early -type population (e.g. 
lAdami et al.1 119981: ICarlberg et all 1 19971) . For example, in the 
Coma cluster, factors up to ~ 1 .4 between the vel ocity disper- 
sion of early- and late-type galaxies are measured (Adami et al. 
U998h . lGotol (120051) suggests that evolved massive galaxies might 
have reduced their velocity by dynamical friction between the 
cluster members through less vigorous tidal interactions (e.g. 
galaxy harassment). Since our sample mainly consists of galax- 
ies with a star-forming SED, the measured velocity dispersion 



might introduce another overestimation of the cluster mass. 
Many approaches to account for the interloper-problem are dis- 
cussed in literature. The above calculation of the virial mass 
(Equation© is an idealised estimate, for the case of perfect iden- 
tification of virialised cluster members. Interloper galaxies lie 
close to the cluster in projected distance and in redshift but are 
in fact members of a different halo or substructure. Restricting 
samples to elliptical galaxies, which unfortunately is not possi- 
ble in o ur case, can redu ce but not alleviate the problem com- 
pletely dWhite et al|2010h . 

Woit ak et al.l (120071) discuss different methods of interloper 
treatment and test them in detail using the results of cos- 
mological N-body simulatio ns.Th ey find that t he m ethods of 
Iden Hartog & Katgertl d!996|) and iPerea et all d 19901) are the 
most efficient ones and achieve an identification rate of gravi- 
tationally unbound galaxies of 60 - 70% at negligib le rates of 
false identifications. Here we adopt the method of IPerea et al.l 
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Table 2. Dynamical parameters for the 4 subclusters using only "virialised"-cluster galaxies. The errors are computed via bootstrap 
simulations (2000 replications). R pv denotes the projected virial radius, M v is the virial mass inside the 3D-virial radius r v = |/? pv . 
p(r„) is the mean interior density inside the virial radius. r2oo is the radius where the mean interior density is 200 p c . M200 is the 
corresp ond ing mass. As a comp arison the corresponding values derived from weak-lensing analysis (H08) following lHevm ans et al. 
d2008l) and lBiviano et alJ d2006l) . 



(1990), which relies on the iterative removal of galaxies whose 
absence in the sample causes the biggest change in the virial and 
projected mass estimator. 

Using only the virialised galaxies identified in this way, we cal- 
culate the dynamical parameters of A901/902 (see Table [2]). 
Also shown are the correspondin g paramete rs derived from the 
weak-lensing analysis dHevmans et alj|2~008l two-halo model). 
Since the authors there define the virial radius r20(w as the ra- 
dius where the mass density of the halo is equal to 200-times 
the critical matter density Q. m (z)p c (z) at the redshift of the clus- 
ter, we first have to extrapolate the mass M200 to our definition 
of f"2oo (Equation |5j. To get an estimation of the velocity disper- 
si on of each subcluster , we then use the M200-O" relation found 
bv lBiviano et al.1 d2006l) . 

Although the applied interloper removal algorithm does lower 
the calculated virial masses, they are partly still significantly 
higher than the values derived from weak-lensing analysis. This 
overestimation is a further indication of the non-virialised state 
of the subclusters and the high significance for substructures ob- 
tained using the DS-test. Hence, the dynamical parameters of the 
dark matter analysis are used hereafter. 

We emphasise that the results of the following section would not 
change significantly if the virial theorem based dynamical pa- 
rameters are used. 

Cluster membership is adopted from Section 3.1, since galax- 
ies residing in the neighbourhood of a cluster centre in terms of 
projected position and velocity can be subject to cluster-specific 
interactions, even though they might be in a non-virialised state. 



4. Kinematic analysis of galaxies 

In this section we analyse the kinematics of our galaxy sample. 
We derive spatially resolved rotation curves from emission lines. 
Since these rotation curves trace the gaseous disk of a galaxy, 
asymmetries in such a curve can be a measure for kinematic 
distortions due to environmental effects like ram-pressure 
stripping, mergers, galaxy harassment or tidal interactions with 
the cluster potent ial. 

iKron berger et aT| (120081) find in combined N- 
body/hydrodynamic simulations that the collision-less stellar 
disk is not affected by ram pressure and only the gaseous 
component shows disturbed kinematics. Asymmetries in the gas 
kinematics combined with a lack of distortions in the stellar 
disk could therefore be a strong indic a tor fo r ram-pressure 
stripping. Furthermore IKronberger et"aD (12008 ) show that the 
effects of ram-pressure stripping depend on the orientation 
of a disk galaxy's plane with respect to its motion through 
the ICM. Characteristic features are a mismatch between the 



kinematic and the luminous centre, declining outer parts of the 
rotation curve or even distortions in the inner parts. However, 
ram pressure was found not to affect the stellar disk of a galaxy 
significantly. 

In this context we want to take a closer l ook at the populati on of 
"dusty red" galaxies first introduced by IWolf et aO (120051) and 
investigate to what extent they differ from "blue cloud" galaxies. 



4.1. Rotation-curve extraction 



First we inspect the one-dimensional spectrum for usable emis- 
sion lines and define a quality parameter ranging from (= no 
detectable emission line) to 3 (strong emission line). Then we 
compute the average in wavelength space of the emission line 
centred 2D-spectrum to get its profile along the spatial axis. Next 
we fit a Gaussian to that profile and derive the luminous centre 
to within 0. 1 arcsec . In cases of galaxies with peculiar light pro- 
files, we redefine the centre manually. The luminous centre thus 
obtained is used to define the kinematic center of a given galaxy. 
Prior to the emission line fitting we average three neighbouring 
rows to enhance the S/N. In the case of very weak emission lines, 
we set this boxcar filter to five rows, corresponding to one arc- 
second (spatial scale = Of.' '205/ pixel). 

We fit a single Gaussian profile to the emission lines 
[O III] A4959A, [O III] J5007A and H/3A486lA and two 
Gaussians to the [O II] /i3726/3729 doublet. For about half of 
the spectra the single profiles to the [O II] doublet yield smaller 
errors and a larger radial extraction range for the rotation curve. 
Subsequently, we apply automatic row-by-row fits to the emis- 
sion line at each position along a galaxy's major axis. The red- 
and blueshifts (with respect to the kinematic centre) along the 
spectral axis can be transformed into observed line-of-sight ro- 
tation velo cities and define an observed rotation curve. For more 
details see Boh m et al.1 d2004l) . 

These rotation curves are not yet corrected for the disk inclina- 
tion 2. However, in this analysis we are only interested in the 
quantitative rotation-curve shape, not the normalisation of the 
rotation velocity. 

We were able to extract rotation curves from at least one emis- 
sion line in 182 different galaxies. Of these, 86 are field galaxies 
and 96 are cluster galaxies. We show an example rotation curve 
with the corresponding 2D emission line spectrum in Figure |5] 
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Fig. 5. Upper panel: The [O II] doublet of an example galaxy. 
Lower panel: Rotation curve extracted from the above image. 
The receding and approaching parts of the rotation curve are 
symmetric. This is an indicator for undistorted kinematics. 



4.2. Rotation-curve asymmetry 

To quantify the kinematic dis tortions of a galaxy we follow the 
approach of D ale et al.l (J2001 ) and use a similar measure for the 
rotation-curve asymmetry. We choose 



A 



Mr,-) + v(-n)\ 
Voffr;) + tf(~n) 



1 y Mr,-)| + H-n)\ 

2 V jo- 2 D (n) + o- 2 v (-n) 



(6) 



The total area between the kinematically folded, approaching 
and receding side is normalised to the average area under the 
rotation curve. Additionally, the contribution of each velocity 
pair (Mr,), v(— r,-)) is weighted by its error (cr„(r,-), <x„(-r,-)). The 
error cr A of the asymmetry is calculated via error propagation of 
the velocity errors o"„(r,-). 

This measure is primarily sensitive to the outer parts of a 
rotation curve, where the faster velocities imply larger absolute 
differences in a velocity pair, while the errors remain the same. 
Accordingly, our asymmetry measure is very sensitive to offsets 
between the kinematic and luminous centre. 
The kinematic centre of a galaxy is the position which by 
definition has no rotational Doppler shift. Since changes in the 
kinematic centre lead to different asymmetry values, we min- 
imise A by allowing shifts of the assumed kinematic centre by 
up to ± 1 .5 pixels with respect to the luminous centre. The centre 
position that results in the smallest asymmetry value is then set 
to be the final kinematic centre for all remaining analysis. Also, 
we inspect the results of the minimisation algorithm visually 
and treat galaxies with peculiar rotation curves or with a small 



radial extent individually. In some cases we reduce the allowed 
shift manually or even set it to 0. 

For undistorted galaxies this minimising algorithm results in 
symmetric outer parts of the rotation curve. Figure |6] shows 
examples of rotation curves with different asymmetry values. 
Seeing correlates the error bars of adjacent points and it can 
smooth out rotation-curve asymmetries, especially in smaller 
galaxies. However, since we do only internal comparisons and 
focus mostly on asymmetries at large galacto-centric radii, the 
impact of these effects on our analysis should be small. 
In Figure [7] we show a histogram o f the offsets between the 
kinem atic and the luminous centres. Men des de Oliveira et al.l 
(l2003h adopt the same technique to determine the offset between 
the kinematic and luminous centre for 25 Hickson compact 
group galaxies. The average distance between the kinematic 
and luminous centre is higher in their compact group galaxies 
(median = 521pc) than in our cluster galaxies (median = 313pc). 
As a further test, we determine a visual asymmetry parameter 




400 600 

offset [pc] 



Fig. 7. Histogram of the offset between the luminous and kine- 
matic centres, determined by minimising the rotation asymme- 
try. 



^visual for each visible emission line, and set it to or 1 
according to whether a rotation curve is visually considered 
undistorted or distorted. Figure |8]shows the correlation between 
the computed asymmetry index A and the visual index A v i sua i, 
which broadly agree. 

If for a given galaxy, rotation curves and subsequently asym- 
metry values A can be determined from more than one emis- 
sion line, we calculate the final asymmetry value via an error- 
weighted average of the values of each emission line. 



4.3. Proxies for environment 

One way to define a proxy for th e cluster environmen t is the 
separation variable s introduced by Carlberg et al1 (ll997l) . which 
uses normalised coordinates and is based on the projected posi- 
tion and velocity of a galaxy relative to the cluster centre. Here, 
we scale velocity differences with respect to each subcluster cen- 
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Fig. 6. Rotation curves with different degrees of kinematic asymmetry. From top/left to bottom/right the asymmetry index A is 
increasing. 



tre to cr sc i and scale projected radial coordinates to r2oo, cl , tne 
radius where the mean interior density is 200p c (z) (see Table|2|: 

r s = — , (7) 

r 200 sd 

frest /Q x 

V s = , (8) 

0"scl 

where r is the projected radius from the subcluster centre and 
y res t is the peculiar rest-frame velocity of a galaxy residing in 



a certain subcluster. Such a scaling allows to compare galaxies 
from different subclusters. The separation s is then given by: 

s 2 = rj + >?s (9) 

For analysing asymmetry relations we combine objects in bins 
and compute for each bin an error-weighted average of the in- 
dividual asymmetry values. This error-weighting is necessary 
since lower quality rotation curves with large error bars would 
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Fig. 8. The visual rotational asymmetry parameter A v i sua i is plot- 
ted vs. the computed asymmetry index A. A v ; sua i is or 1 accord- 
ing to whether a rotation curve is considered as undistorted or 
distorted. The 182 galaxies are combined into bins and than an 
average is calculated (bin size = 45). 



otherwise bias the asymmetry indices to higher values. The er- 
rors of the mean values are calculated taking the maximum of 
the "intrinsic error" (error from the A-values themselves) and 
the "extrinsic error" (error from the scatter of the A-values in 
the bin). The former is usually considered to be the error in the 
mean and is, for most bins, smaller than the error derived from 
the scatter around the mean. So, adopting the latter for the error 
bars is a conservative approach and commonly referred to as a 
"correction for overdispersion". 

4.4. Test of selection and segregation effects 

Before investigating any influence of environment on the 
rotation-curve asymmetry, we try to consider potential conse- 
quences of selection and segregation effects. 
The well-examined morphology-density relation is expected to 
be present in our sample. Hence we should detect a mass and 
luminosity segregation, such that more massive and luminous 
galaxies reside preferentially near the cluster centres. We rely 
on absolut e magnitudes and stellar masses from the STAGES 
catalogue (Grav et al. 2009b) (see Figure|9), and find that stellar 
mass and V-band luminosity increase towards smaller cluster- 
centric radii. Below a radius r s < 0.4 no spiral galaxies are con- 
tained in our sample. This is no surprise since we selected only 
spiral galaxies with a star-forming SED, which are rare in the 
cluster cores. 

Figure [10] shows the asymmetry index A as a function of stel- 
lar mass M* and absolute V-band magnitude V a b s for our sam- 
ple of 96 cluster galaxies. There is a tendency for low-mass 
galaxies to have slightly higher rotation-curve asymmetries. The 
Spearman rank-correlation test, however, shows no significance 
(p = -0.06, p = 41.5%). Low-luminosity galaxies seem to be 
slightly more distorted than galaxies with a high V-band magni- 
tude, but also here we find no significant Spearman-test results 
(p = 0.10, p = 36.2%). Since there are anti-correlations they 
might slightly reduce the average A at small cluster-centric radii. 
However, it is unlikely that these weak trends have a biasing ef- 
fect on our results. 




Fig. 9. Stellar mass M» (blue squares, solid line) and the absolute 
V-band magnitude V a t, s (red circles, dashed line) as a function of 
the cluster-centric radius r s for our sample of 96 cluster galaxies 
(bin size = 32). 
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Fig. 10. The asymmetry index A as a function of stellar mass M« 
(blue squares, solid line) and absolute V-band magnitude V a b s 
(red circles, dashed line) for our sample of 96 cluster galaxies 
(bin size = 32). 



Our original source catalogue comprises ~ 320 cluster galax- 
ies within the selection constraints. Slits could be placed on 103 
cluster members. We could successfully extract rotation curves 
of 96 of these cluster galaxies (success rate ~ 93%, i.e. only 
7 cluster galaxies selected for spectroscopy yield no rotation 
curve). Our sub-sample might not be representative of the whole 
cluster population and differ significantly in some of the prop- 
erties. There could also be regions in the cluster were the spec- 
troscopic completeness is higher or lower. Figure [TT] shows the 
normalised histograms of the cluster-centric distance r s for the 
total source catalogue and for our spectroscopic sub-sample. The 
Kolmogorov-Smirnov statistic, which tests whether two samples 
are drawn from the same distribution, gives a p-value of 47%, i.e. 
we cannot reject the null hypothesis that the distributions are the 
same. Furthermore, the Mann-Whitney-U test, which assesses 
whether one of two groups tends to have larger values than the 
other, is not significant (p-value = 29%). We get similar results 
if we divide our sample into two luminosity bins or stellar mass 
bins. 
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Fig. 11. Normalised histograms of the cluster-centric radius r s 
for the total target sample (blue, 324 galaxies) and our spectro- 
scopically observed VIMOS sample (green, 103 galaxies). 



Additionally, we investigate the luminosity, stellar mass, star- 
formation rate and Sersic index distribution of each sample. In 
every case the Kolmogorov-Smirnov test as well as the Mann- 
Whitney-U test yield non significant differences between source 
catalogue and spectroscopic sample. Also the variation of afore- 
mentioned parameters with scaled cluster-centric distance r s is 
the same within the error bars. 

In summary, these tests clearly indicate that mass (or luminosity) 
segregation, sample selection effects, spectroscopic incomplete- 
ness and other potential biases are unlikely to affect our results 
significantly. 

4.5. Influence of environment 

As outlined in the introduction, cluster environment plays an im- 
portant role in the evolution of galaxies. If cluster environment 
affects the kinematics of galaxies in general or the gaseous disk 
of spirals in particular it should be related somehow to the asym- 
metry of a rotation curve. 

Ram pressure is the force exerted on a galaxy when it moves 
through the hot (~ 10 8 K) and dens e (~ 1CT 3 atoms/cm 3 ) intra- 
cluster medium (lGunn&Gottll972l) . It is stronger at higher rela- 
tive velocities and at smaller cluster-centric radii, where the ICM 
densi ty is higher (proportional to piavi * i^sf e.g. lGunn & Gottl 
1 1 972b . A wide range of numerical simulations and observations 
demonstrate that ram pressure in clusters is sufficient to strip 
part of the gas from galaxies. The stellar disk, on the other hand, 
is not disturbed by the I CM, leaving galax i es with a truncated 



atomic gas (HI) disk ( Boselli & Gavazzi 2006; Chung et al 



2008; Ke nnev et alJ 12004 iKronberge r et all 120081: lOuilis et al 
2000) 



Interactions between a galaxy and the cluster that affect the 
gaseous halo (aka "starvation" or "strangulation") can pre- 
vent the replenishment of the gas disk. A fter consumption of 
the latter, star form ation is quenched (e.g. iBalogh et al. 2000; 
lLarson et al . 1980). In general, such interactions are expected to 
be more efficient in the dense cluster core. 
There are also gravitational galaxy-galaxy interactions, that dis- 
turb the stellar disk of a galaxy significantly. Frequent high- 
speed encounters (harassment) heat the stellar component which 
leads to a more compact morphology due to a decreased rota- 



tional velocity and a thicker bulg e due to an increased veloc- 
ity dispersion (Moore et al In cluster cores mergers are 
rare due to the high relative speeds of galaxie s. Instead, they are 
more likely in the less dense cluster outskirts dBoselli & Gavazzil 
2006), where motions are slower. 

In Figure [12] we show the two asymmetry indices A and A v i sua i 
as a function of the cluster centric radius r s (left panel). We split 
the 96 cluster galaxies in 4 bins with 24 galaxies each and com- 
pute the average asymmetry values in each bin (error-weighted 
in the case of A). The A v i sua i values are rescaled using the rela- 
tion in Figure [8] Both asymmetry indices peak at intermediate 
cluster-centric radii. The Spearman test shows a significant anti- 
correlation (p = -0.42, p = 1.7%). 

In the right panel we plot the asymmetries versus the peculiar 
rest-frame velocity v s . The asymmetry tends to be higher for low 
peculiar rest-frame velocities, although the correlation is not sig- 
nificant (p = -0.22, p = 18.8%), which indicate a mixture of 
different processe s. 

iHeiderman et alJ d2009) find that tidal interactions and mergers 
in the A90 1/902 cluster system dominate in intermediate density 
regions between the cluster core and the cluster outskirts. The 
lack of mergers in the innermost regions is most likely due to the 
large velocity dispersion. Tidal interactions preferably occur at 
low v s values, where the encounters between galaxies last longer 
and the tidal effects can accumulate. In contrast, ram-pressure 
stripping is stronger at higher rest-frame velocities and at smaller 
cluster-cen tric radii, where the ICM-density is higher. 
Vogt et al.l d2004l) investigate optical and H I properties of cluster 
and field spirals, and distinguish four different groups of galax- 
ies. Most field galaxies form a normal population with no pe- 
culiar properties, however, they also find spirals with different 
kinds of peculiarities, which could represent different evolution- 
ary stages in the transformation of field galaxies into cluster SO 
galaxies. Asymmetric spirals have distortions mainly on the lead- 
ing edge of the gas disk. They are probably falling inwards and 
are on their first transit through the cluster. They are assumed 
to have less vigorous interactions with the cluster like e.g. ram- 
pressure stripping. Further classes of stripped (HI deficient; SF 
confined to 3 disk scale lengths) and quenched (H„ in absorp- 
tion) galaxies might have already passed through the core and 
could be transiti onal stages in the transformation process into an 
SO morphology. lVogt et alJ d2004) conclude that gas stripping is 
a significant process in the evolution of spiral galaxies falling 
into a cluster and at least contributes to the morphological trans- 
formation of spiral galaxies into cluster SO galaxies. 
Since a combination of distortions in the rotation curve with a 
lack of asymmetries in t he stellar disk could be an indicator for 
ram-pressure stripping (Kronber ger et al.l 12008 ) we exploit V- 
band HST-images of our sample galaxies to quantify their mor- 
phological asymmetry. 

We here adopt an asymmetry me asure A morp h, which i s fre- 
quently used in the literature (e.g. IConselice et all 12000). It is 
calculated by subtracting a galaxy image 7)8o°, rotated by 180°, 
from the original image /: 



*morph 



Z,-|/-/l80°l 



E; \B - Z?i8o° 



(10) 



B and Biso° are defined analogously to account for the contri- 
bution from sky noise. B and Bigo° are blank sky regions of the 
same size as the original galaxy image. The sums are calculated 
over all pixels within the lcr-isophote of a galaxy. To avoid an 
influence from the determination of the centre of a galaxy, we 
compute a minimisation by shifting the galaxy image on a pixel 
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Fig. 12. The two asymmetry indices A (blue squares, solid line) and A v i sua i (red circles, dashed line) versus cluster centric radius r s 
(left panel) and peculiar rest-frame velocity v s (right panel). The asymmetry is significantly higher at intermediate radii and tends to 
be slightly higher at lower velocities. (96 cluster galaxies; binsize = 24). 
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Fig. 13. The left panel shows the rotation-curve asymmetry A vs. the morphological asymmetry A m0 rph- The black dashed vertical line 
at the median A m0 rph =0.18 divides our sample of 96 cluster galaxies into a morphologically distorted group (MDG) (Amorph > 0. 1 8) 
and a morphologically undistorted group (MUG) (A m0 rph < 0.18) (binsize=16). As a comparison we also show the 86 field galaxies 
(binsize=18). The right panel shows the ratio between the fraction of kinematically distorted cluster /cluster and field / ne id galaxies 
for three different morphological asymmetry bins. 



grid, resulting in t ypical shifts of 2-3 pixels (= 0V06 - 0'. , 09) 
(Boh m et al.1l2012h . 

Figure [13] (left panel) shows the rotation-curve asymmetry A 
vs. the morphological asymmetry A morp h for cluster and field 
galaxies. For several interaction processes (like galaxy-galaxy 
interactions, harassment) we expect coinciding morphological 
and kinematic distortions. However, we find just a slight 
increase in the rotation-curve asymmetry A towards higher 
morphological asymmetry A morp h. On the other hand there 
seems to be a significant fraction of morphologically undis- 
torted but kinematically distorted cluster galaxies. The threshold 
Amorph =0.18 (vertical line in Figure [U} is the median of our 



sample and divides the 96 cluster galaxies into morphologi- 
cally distorted (MDGs) (A morprl > 0.18) and morphologically 
undistorted galaxies (MUGs) (A mor ph < 0.18). Applying a 
Spearman rank-order correlation test to the class of MUGs, we 
find a significant anticorrelation (p = -0.41, p = 6.8%). On the 
other hand we get no significant result for the class of MDGs 
(p = 0.05,/? = 43.8%). 

As a test we apply the same analysis to the 86 field galaxies. 
In contrast to the cluster sample, the peak of A for low A mc . lpn 
is missing. Moreover, the Spearman test above and below the 
median gives no significant result for the field galaxies (MUG: 
p = 0.05, p = 43.6%; MDG: p = 0.24, p = 26.2%). This 
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Fig. 14. The visual rotation-curve asymmetry vs. the morphological asymmetry A morp h for cluster (binsize=16) and field (bin- 
size=17) galaxies. A v i sua i (left panel) classifies galaxy kinematics as undistorted (0) or distorted (1), while A v ; sua i2 (right panel) 
differentiates between three classes: undistorted (0), slightly distorted (1) and heavily distorted (2). 



indicates that the peak at high A and low A morp h for cluster 
galaxies results from a cluster-specific interaction process. 
MDGs might be dominated by tidal interactions while morpho- 
logically undistorted cluster galaxies with perturbed rotation 
curves, might be subject to ram-pressure stripping, an effect that 
is not expected in the field environment. 

Figure [14] (left panel) is similar to Figure Qj] but uses the 
visual rotation-curve asymmetry A v ; sua i. Here, in the cluster the 
fraction of distorted galaxies (63 + 5%) is significantly higher 
than in the field sample (36 ± 5%). Moreover, morphologically 
distorted field galaxies are more likely to also show disturbed 
kinematics (p = 0.25, p = 1.3%). In other words, distortions 
of the gaseous disk in the field are preferentially accompanied 
by disturbances of the stellar disk. This might be due to the 
low relative velocities of galaxies in the field, making tidal 
interactions more efficient. Interestingly, the anticorrelation of 
MUGs in the cluster is no longer significant when considering 
the visual asymmetry parameter. This might suggest that, for 
MUGs, it is not the fraction of kinematically distorted galaxies 
that increases significantly, but rather the intensity of their 
distortions. To verify this conclusion we introduce a refined 
classification parameter, A v ; sua i2 (see Figure[l4] right panel), that 
differentiates between three classes of kinematics: undistorted 
(0), slightly distorted (1) and heavily distorted (2). Using 
this finer scheme we find again a significant anticorrelation 
for morphological undistorted cluster galaxies (p = -0.24, 
p = 4.5%). Figure [13] (right panel) clarifies this conclusion. 
Here we show the ratio between the fraction of kinematically 
distorted cluster galaxies, /duster, and field galaxies, /field, for 
three different morphological asymmetry bins. For morphologi- 
cally undistorted cluster galaxies a mechanism comes into play 
that increases the kinematic asymmetry. 

To summarise this part of our analysis, we found evidence for 
a cluster-specific interaction process that only effects the gas 
kinematics, not the stellar morphology. Ram-pressure stripping 



is, therefore, a prime candidate. 

To illustrate different asymmetry classes, Figure [T5lshows HST 
V-band images with increasing degrees of A morp h and their 
corresponding rotation curves. 

In Figure [16] we show the asymmetry of MUGs and MDGs as 
a function of r s and v s , respectively. We find clear evidence 
for important galaxy populations with distinct characteristics. 
MUGs might be affected by ram-pressure stripping, since they 
exhibit their asymmetry values are anticorrelated with r s for 
r s <> 1.1 (p = -0.85 and p = 0.4%) and positively correlated 
with their rest-frame velocity (p = 0.34, p = 13.0%). Both the 
high ICM density expected at low r s and high relative velocity 
will boost ram pressure and, consequently, produce larger 
kinematic asymmetries. 

On the other hand, MDGs are most likely dominated by 
mergers and tidal interactions. The asymmetry index A peaks at 
intermediate radii (r s ~ 1.0) and is higher for lower rest-frame 
velocities. The Spearman test at greater radii confirms a signif- 
icant anticorrelation (p = -0.45, p = 3.4%). At intermediate 
radii, the galaxy density is already high and the relative galaxy 
velocities are lower than in the cluster core, making tidal 
interactions more efficient. 



4.6. Dusty red galaxies 

In this section we consider dusty red galaxies (selection crite- 
rion in the public catalogue iGrav et all (l2009bl) : sed_type = 2) 
more closely. Twenty-eight of the 96 cluster galaxies in the cat- 
alogue are assigned that SED type. Although they are actively 
star-forming, these galaxies appear red. Their colours are a com- 
bined effect of i) an (on average) four times lower SFR than 
blue cloud galaxies, and ii) intrinsic dust extinction. Having ex- 
tinction levels similar to blue cloud galaxies, dusty red galaxies 
can be understood as the low specific-SFR tail of the blue cloud 
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Fig. 15. Galaxies with different degrees of morphological asymmetry, A morp h, and their corresponding rotation curves. 



dWolfet alJl2003h . While a fit to the colour-magnitude relation 
(CMR) of the red sequence and a parallel cut on its blue side 
(0.25 mag) are used to differentiate between red and blue galax- 
ies, red galaxies with a higher dust red dening then £b-v - 0.1 
are called dusty red (Wolf et al. 2005). They appear predomi- 
nantly in the medium-density outskirts of clusters, while being 
rare in higher and lower density regions. IWolf et all (120051) sug- 
gest two main possible origins. In one scenario, dusty red galax- 
ies representa a transition stage between blue field and cluster 
SO galaxies. The other possibility is that a dusty red galaxy is 
the product of a minor merger between an old red cluster galaxy 
and an infalling bl ue field galaxy. 

Uaffeetal] (12011 ah conclude from studies of the intrinsic colour 
scatter, where they find no significant evolution up to z ~ 
0.8, that most cluster elliptical and SO galaxies have already 
joined the red sequence when the "final" morphology is estab- 



lished. Bamfor det al.l (|2009) detect a considerable fraction of 
spiral galaxies with red colours in denser environments, and a 
substantial population of early-type galaxies with blue colours 
in low-density environments. This implies a dependence of 
colour on environment beyond the morphology-density relation. 
Consequently, the transformation of galaxies from blue to red 
colours must occur on shorter time-scales than any transforma- 
tion from spiral to early-type morphology . These findings favour 
the first scenario proposed bv lWolf et al.l (120051) . in which dusty 
red galaxies are progenitors of SO galaxies. 
We now investigate whether there are differences in the proper- 
ties of dusty red and blue cloud galaxies. The dusty red fraction 
in MUGs is about 8% higher than in MDG, indicating that the 
morphological asymmetry is slightly lower. Similar to FigurefTH 
we analyse the rotation-curve asymmetry A as a function of the 
morphological asymmetry A moip h, for each SED type. Figure [T71 
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Fig. 16. The rotation-curve asymmetry A vs. scaled cluster-centric radius r s (left panel) and scaled rest-frame velocity v s (right panel) 
for MDGs (A moip h > 0.18, red circles and dashed line) and MUGs (A morp h < 0.18, blue squares and solid line). Left panel: MDGs 
have the highest asymmetries at intermediate radii, where tidal interactions are more likely, while MUGs have higher asymmetries 
near the centre. Right panel: the asymmetry of MDGs peaks at low rest-frame velocities, while MUGs have higher asymmetries for 
higher values of v s (binsize=16). 
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Fig. 17. Same as Figure [13] but considering 28 dusty red (red 
circles, dashed line; binsize=9) and 68 blue cloud (blue squares, 
solid line; binzise=12) cluster galaxies, separately. 



shows that there is a dusty red galaxy population with strong 
kinematic distortions but no significant morphological distur- 
bances. The Spearman test for morphological undistorted dusty 
red galaxies (A morp h < 0.18) confirms a significant anticorrela- 
tion (p = -0.68, p = 2.5%), while the correlation for blue cloud 
galaxies is weaker and less significant (p = -0.40, p = 17.6%). 
This favours the scenario that dusty red galaxies are more ex- 
posed to ram-pressure stripping than blue cloud galaxies. We 
stress that many dusty red galaxies with high values of A morp h 
are seen at high disk inclination angles, i.e. they are "pseudo- 
asymmetric" due to the prominent dust disk. To better under- 



stand Figure [T7] we investigate the spatial distribution of the 
emission line flux, which originates from the gaseous disk of the 
galaxy and is therefore much more likely to b e affected by ram- 
press ure stripping than the stellar disk (see e.g. iKronberger et aH 
2008). A measure for the concentration of the line emission is 
the spectroscopic scale length ta^, determined via an exponen- 
tial fit to the emission line flux / along the spatial axis of the 2D 
spectra: 



I(r) = I e~ 



(11) 



Prior to fitting the data, we subtracted the continuum emission 
(estimated red- and blueward of the emission line) to reduce flux 
contamination from the stellar disk. In the wavelength direction 
we averaged flux within a 2 A interval around the emission line. 
The centre of the galaxy was defined to be located at the flux 
maximum. In a first step, we performed an exponential fit on 
each side of this maximum. These values were later used to in- 
fer the gas disk asymmetry. 

However, to derive the total scale length of a galaxy we used a 
more sophisticated fitting process, in which we considered the 
effects of the slit width (in the radio bands often referred to as 
"beam smearing"), the Point Spread Function (PSF), and the in- 
clination, i, of the galaxy with respect to the line of sight. All 
three effects would cause an overestimation of the scale length 
when unaccounted for. 

To correct for the projection effect due to the inclination of the 
disk, the 2D exponential profile was modified in the following 

way: 



I(x, y) = I exp( ) 



I'd 



(12) 



where (xo, yo) are the coordinates of the galaxy center and 
q denotes the ratio of the apparent minor and major axes of 
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Fig. 20. The leading edge ratio /e ra tio (left panel) and scale length ratio r daAo (right panel), respectively, are shown as a function of 
the projected infall-angle <p for blue cloud (blue squares, solid line) and dusty red galaxies (red circles, dashed line). Left panel: 
the leading edge ratio decreases as the infall becomes more edge-on. Right panel: the scale length ratio declines for a more face-on 
infall (binsize: blue clouds=15, dusty reds =10). 




Fig. 18. An illustration of the flux integration perpendicular to 
the slit direction. The exponential profile is projected due to the 
inclination, i, of the disk (here i = 40°). Hence, the isocontours 
have elliptical shapes. The galaxy and slit dimensions shown 
here are typical of our data. In the upper part, the integrated and 
PSF-convolved profiles (solid line) are compared with an expo- 
nential function (dashed line). 



the galaxy. The values of these axial ratio s were previously 
determined using the SExtractor package dBertin & Arnoutsl 
1996) on HST/ACS V-band images, and were adopted from the 
STAGES master catalogue (iGrav et al.ll2009bl) . 
Since the slits were aligned along the apparent major axes 
of the galaxies, the projected exponential profile (Equation 
[T2l) can be integrated perpendicular to the slit direction to get 
the corresponding ID-profile. Figure [18] illustrates how this 
integration was done. The final fitting function was obtained by 
convolving this ID-profile with the PSF. The PSF was assumed 
to be a Gaussian with a FWHM determined from the mean 
DIMM seeing values during the spectroscopy. 
Finally, we correct the best-fit scale length values for the 



position angle 6, which is the tilt angle of the apparent major 
axis with respect to the spatial axis in the 2D spectrum. 
Distorted galaxies have larger fitting errors in their r derec values. 
We flag flux emission profiles that are too peculiar, or for which 
the fit failed, and do not use them in the following analysis. This 
excludes 25 out of 96 cluster galaxies, and 1 1 out of 86 field 
galaxies. 

The photometric scale length is a measure of the concen- 
tration of the continuum emission, which stems from the stellar 
disk. We derived the valu es of the photo metric scale length using 
the GALFIT package (Peng et al. 2002) on the corresponding 
HST/ACS V-band images. In the course of this, we performed 
a bulge/disk decomposition of the galaxies. We assumed a de 
Vaucouleurs profile for the bulge and an exponential profile for 
the disk. 

We normalised the gas disk scale length to the stellar disk scale 
length, by computing the scale length ratio r datia = r dspc Jr dphot . 
When values from more than one emission line were available, 
we computed an error-weighted average. 

We find no significant differences in the mean scale length ratios 
of different SED types. Dusty red galaxies have a mean r datio of 
0.86 ± 0.08, while for blue cloud galaxies we find a value of 
0.99 + 0.07. All cluster galaxies combined show a mean scale 
length ratio of 0.95 ± 0.05, which is about 25% lower than the 
value for the field galaxies (1.27 ± 0.07). This lower average 
scale length ratio of cluster galaxies is most likely the result of 
quenching in progress (removal of the gaseous halo and faster 
consumption of the gas disk) or ram pressure. Note that our 
cluster galaxies are on average 0.27 dex more massive than 
our field galaxies, a difference already evident in the source 
catalogues. We find no significant correlation between mass 
and scale length ratio, which implies a mechanism beyond the 

morphology-density rel ation. 

This fi nding agrees with Bamfo rd. Milvang- Jensen. & Aragon-Salamanc 

(l2007h . who found that cluster galaxies on aver age have a 25% 

lower scale length ratio than field galaxies. Also lMilvang- Jensen! 

(120031) found higher scale length rati os for field galaxie s than 

for cluster galaxies. In contrast to this, Jaffe et all (1201 lbl) found 

scale length ratios of ~ 0.8, independent of environment. 

Our results suggest that the interactions of the gaseous disk 



16 



B. Bosch et al.: Ram pressure and dusty red galaxies - key factors in the evolution of A901/02 







Cluster centre 


Leading edge j 






Galaxy 







Fig. 19. We define the projected infall angle <p as the angle be- 
tween the apparent major axis of the galaxy and the radius vector 
pointing towards the cluster centre. The galaxy edge closer to the 
cluster centre is assumed to be the leading edge. 



of cluster spirals with the ICM, with other galaxies, or with 
the cluster potential result in more centrally concentrated star 
formation. Since the effect is stronger for dusty red galaxies, it 
is likely that they have been exposed to this environment in a 
stronger way, or at least for a longer period. A more centrally 
concentrated star formation might also be important for the 
growth of the bulge component. Besides the depletion of the 
gas, this is a further requirement for the transformation of a 
spiral galaxy into one with an SO morphology. This scenario 
is supported by the bulge-to-total ratios obtained from the 
bulge/disk decompositions using GALFIT. While field galaxies 
(0.086 ± 0.013) and cluster blue cloud galaxies (0.091 ± 0.017) 
have similarly low bulge-to-total ratios, dusty red galaxies have, 
on average, a ratio of (B/T) = 0.232 ± 0.039. 
The position angle (orientation of the apparent major axis) of a 
given galaxy on the sky can be used to introduce a raw estimate 
of the 3-D angle of the galaxy's movement through the ICM. 
We define this projected "infall-angle" <p as the angle between a 
galaxy's apparent major axis and a vector connecting the galaxy 
and the nearest subcluster center (see Figure [T9V We are aware 
that, for a given object, this angle <p might differ strongly from 
the real 3-D movement angle. However, in a statistical sense, it 
should be reliable to use <p to get a rough estimate of the angle 
under which the ICM acts on an infalling galaxy's gas disk. <p 
spans the interval [0°,90°], where 0° corresponds to edge-on 
and 90° corresponds to face-on infall. 

We also define a "leading edge ratio" /e rat i . This ratio compares 
the gas scale length (determined from the 2-D spectra) on one 
side of the disk to that of the other. We compute the le mt io by 
dividing the gas scale length on the side of the disk that (in 
projection) points towards the cluster center (leading edge) by 
the gas scale length on the side of the disk that (in projection) 
points away from the cluster center (trailing edge). It will be 
interesting to see whether the projected infall angle ip has an 
impact on the /e rat ; , which is a proxy for gas disk asymmetry. 
Figure [20] (left panel) shows the leading edge ratio fe ra tio as a 
function of the projected infall angle <p for blue cloud and dusty 
red galaxies. For more edge-on galaxies (smaller <p) the leading 
edge ratio is reduced, implying that ram pressure is successively 
depleting gas from the leading edge. The Spearman test for 



all cluster galaxies confirms a highly significant correlation 
(p = 0.49, p = 0.3%). Furthermore, this clear dependence 
on the infall angle suggests that most galaxies are indeed on 
their first passage through the cluster. Thi s agrees with th e 
aforementioned class of asymmetric spirals (Vogt et al. 2004). 
The authors find that these galaxies have distortions mainly on 
the leading edge of the gas disk, are probably falling into the 
cluster, and are subject to less vigorous interactions (such as 
ram pressure stripping) with the cluster. 

The right panel of Figure [20] shows the scale length ratio r^ 
as a function of the projected infall angle <p. We see a slight 
reduction of for more face-on galaxies (higher <p). However, 
this trend is not significant (p = -0.14, p - 13.2%). Such a 
reduct ion is seen in simulations. For example. ISteinhauser et al.l 
(2012) found in their combined N-body/hydrodynamic simula- 
tions of ram-pressure stripping that the amount of stripped gas 
is larger when the infall orientation is more face-on. 

Figure [21] shows the scaled rest-frame velocity v s as a func- 
tion of the scaled cluster-centric radius r s . While blue cloud 
galaxies show no significant correlation (p = 0.16, p = 11.2%), 
dusty red galaxies have higher velocities at large radii (p = 0.46, 
p - 0.7%). In contrast to our dusty red and blue cloud sample, 
early type galaxies have, in general, higher velocities at smaller 
radii, due to the deeper gravitational potential well towards the 
centre. This discrepancy might be partly due to a high fraction 
of infalling and "non-virialised" galaxies. The strongly peculiar 
relation for our dusty red galaxies could be explained within the 
framework of ram pressure stripping. If we assume that most of 
the dusty red galaxies are falling into the cluster (which is sup- 
ported by Figure l20l this finding could best be explained as a 
selec tion effect, since ram pressure is proportional to picm * u 2 el 
(e.g. iGunn & Gottlll972l) . Under the assumption that dusty red 
galaxies are undergoing (and/or produced by) strong ram pres- 
sure, galaxies residing in the outskirts of the cluster need a high 
relative velocity to compensate for the low ICM density. On the 
other hand, for dusty red galaxies with high rest-frame velocities 
residing close to the cluster centre, ram pressure might be so high 
that these objects are already quenched, with no detectable star 
formation. Since we have selected our galaxies by a star-forming 
SED, these quenched objects with high velocities at low radii are 
missing in our sample. 

For the Virgo cluster, Chun g et al.l d2008) found stripped H I 
disks outside the core region (> IMpc), some of them show- 
ing signatures of ongoing ram-pressure stripping and provid- 
ing evidence that ram pressure is also effective at such dis- 
tances. Exploiti ng the dark matter pro perties of the weak lensing 
measu rements dHevmans et al.ll2008l) . we follow iMakino et akl 
(119981) and estimate the distribution of the gas in each subcluster. 
Assuming an isothermal /3-model f or the intra-cluster medium 
(ICavaliere & Fusco-Femianolll976l) . we get ICM-density esti- 
mates that are p gas ~ 10" 28 gctrT 3 at 1.5 * r^oo, high enough 
to produce significant ram pressure if galaxies have relative 
velocities around 1000 kms" 1 a s has been sho wn with sim- 
ulation s ISteinhaus er et al.l 120121) . Furthermore, iKennev et al.l 
(2004) suggested that a merging of subclusters (which is the case 
for A90 1/902) can cause large bulk motions and local density 
boosts, raising th e ICM pressure t emp orarily. This is also seen 
in simulations by iTonnesen et al.1 d2007l) . and it facilitates ram- 
pressure stripping in the outskirts. These ICM density boosts 
might not show up on shallow X-ray maps like our XMM data 
ofA901/902. 

At a given relative velocity, the time scale for a galaxy to lose all 
of its gas due to ram pressure will depend on the ICM density. 
That is, it is less likely to observe a relatively fast galaxy that 
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Fig. 21. Scaled peculiar rest-frame velocity v s is shown as a func- 
tion of the scaled cluster-centric radius r s for blue cloud (blue 
squares, solid line) and dusty red galaxies (red circles, dashed 
line). High-velocity dusty red galaxies are predominantly re- 
siding in the cluster outskirts (binsize: blue clouds =22, dusty 
reds=9). The black dotted lines indicate lines of constant ram 
pressure for each subcluster, normalised to v s =r s =l. 



still has some gas disk reservoir closer to the cluster core. This 
could explain the decrease in high velocity dusty red galaxies 
at smaller cluster-centric radii. For blue cloud galaxies this ef- 
fect might be less severe, due to higher gas mass fractions when 
falling into the cluster system. To strengthen this conjecture, the 
black dotted lines in Figure [21] indicate the regime of constant 
ram pressure in the r s -v s plane. The constant is set to be the ram- 
pressure value at r s =v s =l, i.e. const = p gas ('".!=l) * v K i(y s =l) 2 . 
Since the /3-profile scales with 1/r 2 for radii well outside the 
core-radius, u 2 /r 2 = const, or v = const * r. Thus, iso-ram- 
pressure contours are just straight lines in this diagram. 



5. Summary and conclusions 

We presented spectroscopic observations of 96 disk galaxies in 
the Abell 901/902 multiple cluster system, which is located at a 
redshift of z ~ 0.165. Additionally, we compiled a comparison 
sample of 86 field galaxies at similar redshifts. The data were 
taken with the VLT instrument VIMOS. We carried out a 
redshift analysis and estimated dynamical parameters of the 
four subclusters, including their virial masses, virial radii, 
and velocity dispersio ns. For this we used biweight statis 



ana velocity dispersions, tor tms we used biweignt statis- 
tics dBeers et al.l[l 990) an d an i nterloper removal procedure 



suggested by Pereaetal. (1990). We ap plied the Dressler- 
Shectman test (Dre ssier & Shectmanl H988) to identify possible 
substructures in the cluster system. Furthermore, we extracted 
rotation curves from the spatially resolved emission lines and 
analysed distortions in the gaseous disk of the galaxies as well 
as HST/ACS images to quantify peculiarities of the stellar disk. 
We paid particular attention to dusty red galaxies as a possible 
intermediate stage in the transformation of field spirals to cluster 
SO galaxies. The results of this work can be summarised as 
follows: 



and in a non-virialised state. Accounting for interlopers, we 
derived the dynamical properties of the four subclusters. 

2. The fraction of kinematically distorted galaxies is 75% 
higher in the cluster than in the field environment. This 
difference mainly stems from morphologically undistorted 
galaxies, indicating a cluster-specific interaction process 
which only affects the gas kinematics, not the stellar mor- 
phology. Ram-pressure stripping is therefore a prime candi- 
date. 

3. We find two important galaxy populations: 

- Morphologically undistorted galaxies, which show high 
rotation-curve asymmetries at high rest-frame velocities 
and low cluster-centric radii, indicating that this group 
is strongly affected by ram-pressure stripping due to the 
interaction with the intracluster medium. 

- Morphologically distorted galaxies, which are probably 
subject to tidal interactions. They show higher rotation- 
curve asymmetries at intermediate cluster-centric radii, 
and low rest-frame peculiar velocities. 

4. Among the morphologically undistorted galaxies, the dusty 
red galaxies have particularly high rotation-curve asymme- 
tries, suggesting an enhanced effect of ram pressure. 

5. The ratio between gas and stellar scale length is reduced for 
cluster galaxies compared to the field. This indicates that 
disks in the cluster environment lose part of their gas reser- 
voir. 

6. The projected disk orientation of the cluster galaxies seems 
to be correlated with the gas scale length in the sense that: 

- The ratio between the gas and stellar scale length steadily 
increases, going from face-on to edge-on infall. 

- For more edge-on infall, the gas scale length of the lead- 
ing edge is reduced compared to the gas scale length of 
the trailing edge of a galaxy. 

Both findings suggest an influence of ram pressure, and are 
in agreement with numerical simulations. Furthermore, they 
indicate that the spiral galaxies in our sample are in the pro- 
cess of falling into, and on their first transit through, the clus- 
ter. 

7. Dusty red galaxies on average have bulge-to-total ratios a 
factor of 2.5 higher than cluster blue cloud and field disk 
galaxies. This is further evidence that dusty red galaxies rep- 
resent an intermediate stage in the transformation of field 
spirals to cluster lenticulars. 

8. Dusty red galaxies residing at the outskirts of the cluster have 
a higher rest-frame velocity. This might be indirect evidence 
for strong ram pressure producing this SED-type, as high ve- 
locities compensate the low ICM density at these radii. 

9. Galaxies with high velocities near the cluster centre are prob- 
ably already quenched and are scarce in our sample due to 
our selection for galaxies with a star-forming SED. 

We interpret the aforementioned results as indications that ram- 
pressure stripping is a key factor for galaxy transformations in 
the multiple cluster system A90 1/902, and that dusty red galax- 
ies could be a transition state from spiral to SO galaxies. 
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